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We investigate quantum phase transitions (QPTs) in spin chain systems characterized by local Hamiltonians 
with matrix product ground states. We show how to theoretically engineer such QPT points between states with 
predetermined properties. While some of the characteristics of these transitions are familiar, like the appearance 
of singularities in the thermodynamic limit, diverging correlation length, and vanishing energy gap, others differ 
from the standard paradigm: In particular, the ground state energy remains analytic, and the entanglement 
entropy of a half-chain stays finite. Examples demonstrate that these kinds of transitions can occur at the triple 
point of 'conventional' QPTs. 
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A considerable part of modern condensed matter physics is 
devoted to the study of matter near zero temperature. In partic- 
ular, zero- temperature quantum phase transitions (QPTs) 0], 
as observed in cuprate high-temperature superconductors and 
heavy fermion materials, have attracted enormous attention. 
Although an adaptation of the classical Landau-Ginzburg the- 
ory successfully describes some of these phase transitions, it 
is manifest that this concept is in general too narrow to cover 
all the fascinating aspects of the quantum world 0] ■ A com- 
plete and rigorous quantum mechanical description is, how- 
ever, burdened by the notorious complexity of quantum corre- 
lations in highly entangled many-body systems. 

The fields of condensed matter and quantum information 
theory study the behavior of quantum many-body systems by 
using complementary methodologies. Whereas the typical 
starting point in condensed matter theory is a Hamiltonian, 
from which states emerge as ground states (GSs) or excita- 
tions, quantum information theory deals primarily with quan- 
tum states, from which corresponding Hamiltonians may be 
constructed. For spin chains this point of view can be traced 
back to the seminal works on the AKLT model |3| and finitely 
correlated states fl], and it has recently been successfully re- 
sumed in various works on matrix product states (MPS) @|. 
This led to new powerful numerical algorithms jg, 0] ac- 
companied by a better understanding of their efficiency @L 
and new insights in renormalization group transformations 
fioll and sequential quantum generators full . Clearly, a fruit- 
ful crossfertilization is emerging between these two fields. 

The present work investigates QPTs in systems represented 
by MPS by following the quantum information approach. It 
thus generalizes the findings of which already indi- 

cated the possibility of such transitions in MPS systems, and 
it enables us to theoretically design QPT's in quasi-exactly 
solvable models. In this way, we may engineer QPT points 
between phases whose correlations or symmetries we choose 
a priori. The corresponding orders can be of local type and/or 
of a more subtle hidden non-local character. The main obser- 
vation behind is that, for the systems under consideration, a 
singularity in a (T> — l)-dimensional transfer operator leads 
to a QPT in the corresponding 2?-dimensional quantum sys- 
tem. Although these findings hold for arbitrary T>, we focus 



on MPS in T> = 1, where a general discussion is possible on 
full analytic grounds. 

In fact, every state, in particular every GS, of a finite system 
can be represented as a MPS 01 ■ The power of this rep- 
resentation — and with it the power of Density Matrix Renor- 
malization Group (DMRG) — stems from the fact that in many 
cases a low-dimensional MPS already leads to a very good ap- 
proximation of the state \9]. From such a low-dimensional 
MPS one can in turn construct a parent Hamiltonian from 
which it arises as an exact GS. We will study the dependence 
of correlation functions of such systems on a smoothly vary- 
ing parameter g and show that singularities can appear, which 
are reminiscent to those arising in known examples of QPTs. 
They appear only in the thermodynamic limit and are accom- 
panied by diverging correlation lengths and vanishing energy 
gaps. 

Some of the derived properties do, however, hardly fit 
within the conventional picture of QPTs in V = 1 spin sys- 
tems: First, at the QPT point g = g c the GS energy density 
eo is analytic (Fig. 1). Typically, a non-analyticity in eo(g) is 
used as the defining property of a QPT [l|] . We think, however, 
that the use of the term QPT is well justified by the presence of 
a non-analyticity of any observable quantity (e.g., two-point 
correlations). The second non-typical feature is the fact that 
the entropy of a half-chain remains finite as g — > g c , which re- 
flects the fact that MPS-QPTs cannot be described in terms of 
conformal field theory. In fact, a lot of attention has recently 
been devoted to the entanglement entropy Jl3lfl4ll . resulting 
in the observation that the crossing of a QPT point typically 
coincides with the divergence of this entropy. The discussion 
below shows, however, that this is not the case for MPS QPTs 
in T> = 1. A third peculiarity is, that although the breaking of 
a discrete symmetry can be engineered, MPS QPTs can occur 
without spontaneous symmetry breaking since for g ^ g c the 
GS is typically unique. 

We will start by discussing the relevant properties of MPS 
and show how singularities can appear in the thermodynamic 
limit (size of the chain N — > oo). Then, following the idea 
of |4J we show how the corresponding Hamiltonians can be 
constructed and discuss two examples, with local and non- 
local order, in greater detail. For one of these we show that the 
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MPS-QPT occurs exactly at the triple point of 'conventional' 
QPTs. 

Singularities in MPS — Consider a MPS which is up to 
normalization given by 



J2 tr[A 4l ...A 



(i) 



ti...ijv=l 



where {Aj} is a set of d D x £) matrices, d is the Hilbert 
space dimension corresponding to one site in the chain and D 
is the dimension of the bonds, when we think of the state in 
the valence bond picture. The state in Eq. Q is translational 
invariant on a ring of length N, it has reflection symmetry if 
A = Af, permutation symmetry if [A, Aj] = and time 
reversal symmetry if the Aj's are real. Also other local sym- 
metries (e.g., SU(2), 1i) can be enforced by imposing appro- 
priate constraints on the A 

In the following we will consider systems where the matri- 
ces Ai depend on a single real parameter g. It is important to 
note that if the Ai depend on g in an analytic or continuous 
way, then so will its parent Hamiltonian which we construct 
below. Clearly, we could consider more general cases with 
several parameters gi, g%, . . . . 

Correlation functions for m consecutive sites are given by 



(Si . . . S m ) 



tr[E»- m E Sl ...E s J 
tr[<] 



Es = ^(WMj®^ 



, with (2) 

(3) 



Here Si is any observable acting on the i th site and the bar 
denotes complex conjugation. For simplicity we will focus 
on the generic case where the transfer operator E\ is diag- 
onalizable and non-degenerate for g ^ g c . Taking the ther- 
modynamic limit (N — > oo) only the right \r) and left \r) 
eigenvectors of Ei's largest eigenvalue vi survive in Eq. Q. 
With the normalization (l\r) = 1 this leads to 



(Si...S m ) = (l\E Sl ...E s Jr)/v? 



(4) 



Hence, if for some g — g c there is a level crossing in the 
largest eigenvalues of E%, then there will typically be a dis- 
continuity in the correlation functions (or their derivatives), 
even though the Aj's and with them the E's are analytic in g. 
Needless to say that the same argumentation holds for every 
observable with finite support. 

A trivial example showing that discontinuities of any order 
n are possible is given by (D = d = 2) 



Ai = 



1 
1 + 9 



A 2 = 



g n 




(5) 



Here, all derivatives dg (SjSj+i) of order k < n will be con- 
tinuous at g = g c = 0, whereas the n th -order derivative e.g. 
of (a x ) turns out to be discontinuous. 



Let us now discuss the properties of a general MPS in the 
vicinity of a transition point g c . The decay of two-point cor- 
relations (SiSi+i) can be obtained from Eq. (0 by setting 
Es 2 = ■ ■ ■ = Es m _ 1 = Ei with m = I + 1 and exploiting the 
Jordan decomposition of the transfer operator. This leads to 



(SiS i+l ) - (Si)(S i+l ) 



l-l 



(6) 



where v 2 is the second largest eigenvalue of the transfer opera- 
tor. As the coupling strength approaches its QPT point value, 
g — > g c , we get \v?\ — > \vi\- Then, the correlation length 
£ = 1/ log | i^i /V2 1 diverges and one obtains long-range cor- 
relations at the transition point. Note that despite the diverging 
correlation length there is no power-law decay at the transition 
point (which can be different for V > 1 Q(J). 

Since a lot of attention has recently been devoted to the re- 
lation between criticality and the scaling of the entanglement 
entropy fl3ll we give now an explicit formula for the latter. In 
fact, in recent works QPT points seem to be intimately con- 
nected with a logarithmic diverging behavior of the entropy 
of a block of consecutive spins when considered as a function 
of the block size. However, this is not the case for the class 
of QPT's discussed in the present work. The entropy of an 
asymptotically large block of a MPS can be calculated exactly 
by exploiting the freedom in the A's in order to fix the gauge 
AiA\ = 1, Y^i AlgAi = g, where g is a density matrix 
acting on C D . By the renormalization group arguments of 
M il (III the spectrum of a large block converges to the spectrum 
of q® 2 such that the entropy becomes 2S(g) = — 2trglog 2 g. 
In particular, the entropy never exceeds 2 log D, irrespective 
of how close the system is to a QPT point. This immediately 
implies that states exhibiting a diverging growth of the en- 
tanglement entropy cannot be described by MPS with finite 
D (explaining the difficulties of DMRG for critical systems). 
Moreover, following 11711 the absence of a logarithmic diver- 
gence implies that MPS-QPTs cannot be described in terms of 
conformal field theory (as already indicated by the absence of 
power-law decaying correlations). 

Higher dimensions — Although our main focus is on one- 
dimensional quantum systems {T> = 1) we briefly discuss in 
this interlude how to extend the presented ideas to higher di- 
mensions. 

Note that in the above T> = 1 case the QPT is traced back 
to a level crossing in the operator E 1 , which acts only on a 
single site, i.e., in dimension T> — 1. However, if the trans- 
fer operator itself has a spatial substructure, we are led to a 
QPT in a higher dimensional system. To be more specific 
consider a P-dimensional cubic lattice of size N\ X . . . X N-p 
with periodic boundary conditions. A generalization of ©, 
the so-called projected entangled pair state (PEPS fiH ). is 
then obtained by replacing the bi-linear forms Aj by tensors 
of order 22?, i.e., A : (C D ) mv -> C and the matrix product 
by a tensor contraction according to the edges of the lattice. 
As we can interpret the 2?-dimensional lattice as a chain of 
2? — 1 dimensional systems we can introduce a transfer opera- 
tor Ei for this chain by contracting all N% operators E t 
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FIG. 1: Energies for the ground and first excited states as a func- 
tion of g. In contrast to other QPTs the ground state energy of 
a MPS Hamiltonian remains (by construction) analytic at the QPT 
point g = g c . Nevertheless the spectral gap vanishes and the cor- 
relation length diverges. Moreover, the non-analytic change of the 
ground state at g c is reflected by an observable non-analyticity of 
certain local expectation values. 

on a V — 1 dimensional sub-lattice. In this way we are back 
to the one-dimensional scenario described above. That is, if 
we take N\ — * oo then a level-crossing in the largest eigen- 
value of E' x can give rise to a QPT. Note that this way of con- 
structing QPT points resembles the classical transfer matrix 
method discussed for instance in fiUl . The problem of calcu- 
lating the largest eigenvalue of E[ can numerically be tackled 
by DMRG (PEPS) algorithms for V = 2 (V > 3V Particu- 
lar instances of such transitions were discussed in 1 20], where 
power-law decaying correlations could be determined using 
Monte-Carlo methods. Note that by construction all critical 
ground states obtained in this way obey an area law for the 
entanglement entropy (3- A method for deriving analytic 
results for particular higher dimensional instances will be pro- 
vided elsewhere flrjll . 

The Hamiltonians — Following the works on the AKLT 
model and finitely correlated states one can always construct 
a local Hamiltonian such that a given MPS is its GS: since the 
reduced state density operator pk corresponding to k sites of 
a MPS has at most rank D 2 , it has a nullspace whenever k > 
log D 2 / log d. Therefore, \ip) is the GS of any Hamiltonian 
which is a sum of (local) positive operators supported in that 
nullspace. In particular, it is the GS of the Hamiltonian 

H = Ys T *(Pk), (7) 

i 

with being the projector onto the nullspace of pf. and Tj its 
translation to site i. By construction the GS energy is always 
zero, i.e., it is evidently analytic in g. 

Let us now see in which cases H depends analytically on g 
and, moreover, is such that \ip) is its unique GS for g ^ g c . 
Consider to this end the operator 

R = A® k (t D ®Lo® (k - 1] ®l D )A^ k , (8) 

where w = Ey=i Rfe'l wdA\a,j3) = EJ^k/3K)- 
Note that, if Ai depends analytically on g, R also depends ana- 



lytically on g. It is evident from the valence-bond construction 
of the MPS that in general, range(pk) C range(R). How- 
ever, if both eigenvectors |r) and |Z) have full Schmidt rank, 
then a straightforward calculation shows that range(pk) = 
range(R). Hence, if this Schmidt rank condition is satisfied 
on both sides of g c (which is generically the case), then H (g) 
indeed depends smoothly on g. 

The uniqueness of the GS of Hamiltonians of the form @ 
was discussed in where it was proven that it is unique 

if the largest eigenvalue of E 1 is non-degenerate (i.e., g ^ g c ), 
rank(pk) = D 2 and 

range(pk) ® C d fl C d ® range(pk) = range(pk+i) ■ (9) 

It is also shown there that the latter condition is always satis- 
fied if we replace fc by k + 1, i.e., take H = J^i T i{Pk+i) ■ 

The analyticity of H together with the uniqueness of its GS 
for g 7^ g c immediately imply that a non-analyticity in the 
expectation values can only be caused by a vanishing energy 
gap at g c . 

Note that a degeneracy in the GS is equivalent to a degener- 
acy in Ej . This means that there is no spontaneous symmetry 
breaking in one of the phases, unless \v\\ = for an en- 
tire interval (e.g., for g > g c ). However, for a degenerate E\ 
arbitrary broken discrete symmetries are possible 11511 . 

Examples — After having discussed the general properties 
of MPS QPTs, we will consider some examples in more detail. 
Since the MPS have a very efficient parametrization in terms 
of the matrices Ai we can, by imposing constraints on these 
matrices, engineer systems having desired properties (sym- 
metries, orders, discrete symmetry breaking, etc.). 

1. Three-body interactions: We start by considering the 
case D = d = 2, i.e., A\,Ai being two-by-two matrices. 
By the arguments above every such state has a parent Hamil- 
tonian with local three-body interactions. In fact, many of 
these Hamiltonians are similar to those appearing in (triangu- 
lar) optical lattices l2liE^l . We construct an example with Z2 
symmetry by imposing the existence of a similarity transfor- 
mation which interchanges A\ and A2, i.e., X" 1 A\X = A2 
and X~~ X A2X = A±, This is indeed the case if we choose 

*-(;;)• *-(is)'- ,x -(!s)- (w) 

The corresponding transfer operator E 1 has largest eigenval- 
ues 1 ± g leading to a singularity at g = g c = 0. A straightfor- 
ward calculation shows a discontinuity in the first derivative 
of (afaf +1 )(a = x, y, z), whereas all two-point expectation 
values are continuous. The magnetization in the ^-direction 
can serve as an order-parameter since (cr x ) = 4g/(l + g) 2 
for g > 0, whereas it vanishes for g < 0. 

At the QPT point g c the state is a GHZ state. For g = — 1 
(disordered phase) it is equal to the cluster state I23I1 . and for 
g = 1 (ordered phase) all spins point in the a>direction. 
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The parent Hamiltonian is (by construction) Z2 symmetric: 

H = Y, 2 (-9 2 - 1) °?°f+i - (! + 5) 2 dD 

+ (<?-l)V^i<+2- (12) 

Hence, it is a combination of an Ising interaction with trans- 
verse magnetic field Jilt , and a cluster state Hamiltonian ill 2b . 
Since a constant term was omitted, the GS energy density is 
eo = —2(1 + g 2 ), and one can readily check condition l|9} 
implying that the GS is indeed unique for g 7^ g c = 0. 

The above Hamiltonian can be embedded into a two- 
parameter family 

H{j,h) = -- 2^^- ? — 

+h cr-a- +1 , 

which can be mapped onto a system of non-interacting 
Fermions by a standard Jordan-Wigner transformation. More- 
over, H(j,h) can be mapped onto the XY-model via a duality 
transformation I24IE5I1 . This in turn exhibits 2 nd order QPTs 
on the lines h = ±1 and 7 = for h G (—1,1). The path pa- 
rameterized by g is given by 7 2 + h 2 = 1 (the disorder-line in 
the XY-model). Hence, the MPS transition occurs at the triple 
point of 'conventional' QPTs exhibiting algebraically decay- 
ing correlations and diverging entanglement entropies 1E4I1 . 

2. Two-body interactions: The previous example corre- 
sponded to the case of a local order parameter. Let us fi- 
nally discuss an example with non-local string order. To this 
end, consider D = 2, d = 3, i.e., states of a spin one chain 
which are GSs of nearest-neighbor interactions. The most 
popular MPS in this class is certainly the GS of the spin- 
1 AKLT model, which exhibits a hidden (string) order. In 
fact, this state can be embedded into a one-parameter family 
with MPS QPT If we choose {A,} = { - a z ,a~,ga+} , 
then the transfer operator E\ has eigenvalues —1 and 1 ± g 2 
leading to a diverging correlation length for g — > g c = 0. 
Moreover, the first derivative of (S z ) has a discontinuity at 
g c = 0, where (S z ) — ► 1, i.e., the state becomes ferromag- 
netic. For g = 2 we get the AKLT state and for g = — 2 
a state, which is equivalent to the latter up to a local unitary 
transformation. For g — > ±00 the GS becomes the Neel GHZ 
state (| TIT •••> + ! iU •••>)/\/2. As already shown in 
0] (with a different parametrization and for g < 0) the cor- 
responding Hamiltonian is rotationally symmetric in the XY- 
plane, gapped with a non-degenerate GS (unless g = g c ) and 
has the form 

H = ^(2 + 3 2 )^S ? t+1 +2(S ? t S ? l+1 ) 2 (13) 

i 

+2(4-.g 2 )(5f) 2 + (.g + 2) 2 (^^ +1 ) 2 
+g(g + 2){S!S* +1 ,S i S i+1 } + . 

Note that since E\ (g) = E% (—g) the GSs corresponding to 
±g merely differ by local unitaries. 



Conclusion — Matrix product states provide the perfect 
playground for investigating novel types of quantum phase 
transitions that do not fit in the traditional framework. We 
provided an example of such a QPT at the triple point of con- 
vential phase transitions, and it is easy to construct many more 
such transitions: given two predetermined MPS and associ- 
ated Hamiltonians for which they are the unique ground states, 
we can always construct a one-parameter family of MPS in- 
terpolating between the two. If on the chosen path the trans- 
fer operator exhibits a level crossing in the largest eigenvalue, 
then the system undergoes a QPT. 
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